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Abstract 

We consider two types of non-autonomous 2-periodic Gumovski-Mira difference 
equations. We show that while the corresponding autonomous recurrences are conju- 
gated, the behavior of the sequences generated by the 2-periodic ones differ dramatically: 
in one case the behavior of the sequences is simple (integrable) and in the other case it 
is much more complicated (chaotic). We also present a global study of the integrable 
case that includes which periods appear for the recurrence. 

Keywords: Integrable and chaotic difference equations and maps; rational difference equa- 
tions with periodic coefficients; perturbed twist maps. 
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1 Introduction 



Autonomous difference equations are a classical tool for the modeling of ecological systems, 
see for instance \W\ [TE[ 1^ . One of the modifications applied in the models in order to adapt 
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them to more realistic situations consists in converting the recurrences into non-autonomous 
ones changing one of the constant parameters by a periodic cycle (see [9l [IHl HH EHl HU [2T] 
and references therein). In this situation it is said that the model takes into account 
seasonality. For instance a parameter taking values in a cycle of period 4 could model an 
ecological situation that has different features during spring, summer, autumn or winter. 

In this paper we will consider a very simple autonomous recurrence, which depends on 
a unique parameter and another one which is conjugated to this one. Then we will show 
that changing the constant parameter by a two cycle has a completly different effect in both 
cases. Indeed in one of them the recurrence has an invariant while in the second one the 
recurrence looks as a chaotic one. This phenomenon is very surprising and shows that this 
procedure used in modelling of changing the constant parameters by periodic cycles is quite 
delicate. 

More specifically, we take 



Xn+2 = -Xn + with Ct > 0, (1) 

which is one of the recurrences considered by Gumovski and Mira in [T7]. Performing the 
change of variables x„ = \/a it writes as 

yn+2 = -yn + r, ^"'^^ — wit'^ /? = l/a > 0. (2) 

We will see that when considering the corresponding 2-periodic recurrences 

Xn+2 = -Xn + 2 ' 

and 

I Xn+l 

Xn+2 = -^n+ 2 ' (4) 

Pn -t- -^n+l 

where {an}n is a 2-periodic cycle of positive values and /3„ = l/a„,n > 0, their behaviors 
are completely different. For the second case it is very simple, indeed integrable (that is, it 
has a two-periodic invariant), see Theorem [sj On the other hand, for the first recurrence 
and suitable initial conditions the points {xn}n seem to fill densely many disjoint closed 
intervals. In fact its unfolding (x„, x^+i) presents all the complexity of the perturbed twist 
maps, see Section [3] and also Figure 1. 

More specifically, in Subsection 2.2, we characterize all the sequences generated by the 
recurrence (|4]). They are: 

(i) Constant sequences or periodic sequences with period 2q for some g G N; 

(ii) Sequences such that its adherence is formed by the same orbit plus an accu- 
mulation point; 
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(iii) Sequences {x^jneN such that its adherence is formed by one closed interval; 

(iv) Sequences {x^jneN such that its adherence is formed by two closed intervals. 

With respect to recurrence ([S]) , we will see in Section [s] that there appear the same four 
types of sequences given above plus other ones: 

(v) Sequences such that its adherence is formed by the same orbit plus some 
accumulation points; 

(vi) Sequences {xn}neN such that its adherence is formed by several closed intervals. 

The differences are perhaps more clear when considering the composition maps ([6, 7J) 
associated to the non- autonomous difference equations, which in this case are given by 
Fa2,ai ■= Fa2° Fai and Gp^^p-^ ■■=Gp^o Gp-,, where 

Fa{x,y) = ( j/, -X + ) and Gi3{x,y) = ( y, -x + ^ 



are the corresponding maps associated to the recurrences in ([T]) and ([2]). Notice that for 
instance 

(Xi,X2) > (X2,a;3) > (X3,X4) > (X4,X5) > (X5,X6) > ■ ■ ■ 

where {x^j^gN is the sequence generated by Q and similarly for the recurrence ([3]). 

For all a and /3 the maps Fa and Gp are integrable diffeomorphisms from to M^, 
which preserve area, and the dynamics on the level sets of their corresponding first integrals 
is translation-like ( O [5] ) . In a few words this means that when the level sets are circles it is 
conjugated to a rotation and when they are diffeomorphic to the real line, it is conjugated 
to a translation on this line. 

The curious phenomenon that we present is that while the maps Fa and Gi/a are 
conjugated via the linear change ^(x, y) = {^/ax, \/ay) (that is Gi/a = ° Fa^"^), and 
both maps posses a first integral, the dynamics of the compositions maps are very different: 
while Fa2,ai numerically exhibits all the features of a non- integrable perturbed twist map, 
the map 6*1/0,2,1/01 has a first integral given by 

N 1 2 , 1 2 , 2 2 
V[x,y) = — X H y +x y — xy, 

and its dynamics is translation-like (see Theorem [2]) . These two different behaviors are 
illustrated in Figure 1 below. 
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2 The integrable case 

Consider the non-autonomous recurrence 



Xn^2 = -Xn + r. , 2 ' 

Pn + 

with 

{a > for n = 2k, 
6 > for n = 2k + l, keN, 

The local analysis of the constant recurrences generated by it and the level sets of its 
associated non-autonomous invariant I{x,y,n) = PnX^ + /3„+iy^ — xy + x^y^, has been 
already done in |8] (see also [3] and [H!). In this section we complete the description of the 
global dynamics of the difference equation ([5| and its associated composition map G^^a- 

2.1 Dynamics of the discrete dynamical system generated by G^ a 

The above invariant for ([s]), gives the first integral 

Vb,a(2;, y) ■= y, 2) = ax'^ + hy'^ -xy + x^y^, 

for the the map Gb^a, 

n ( \ f , y {a + y^){y-x{a + y^)) \ 

Gb,a{^,y) = l-X^ ■ 2'~y+'u , 2\2"T7 ? , 2V^ • 

V a + y^ b{a + y^Y + [y - x[a + y^)Y J 

In [8] a full description of the level sets of V^^a is given, they are also described in Theorem 
[2| bellow. In a few words, the level sets are like (a) in Figure 1 when ah > 1/4 and like 
Figure 2 when ah < 1/4. 
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Figure 2: Some orbits of the maps Gh^a with a = 1/4, 6 



2/9. 



To describe the dynamics on each level set we need to apply the result bellow, which 
follows from [S]. 

Proposition 1. Let (f) : U ^ U, U C be a smooth planar area preserving map having 
finitely many fixed points and a smooth first integral V : U ^ M. Then: 

(a) When a connected component of a level set ofV does not contain fixed points of ^, is 
invariant and 

(i) it is diffeomorphic to a circle then <^>, on it, is conjugated to a rotation, 
(a) it is diffeomorphic to the real line then ^, on it, is conjugated to a translation. 

(h) When a connected component of a level set contains fixed points of ^. Then on each 
invariant connected component of this level set minus the fixed points, ^ is conjugated 
to a translation. 

Taking into account the above proposition and the topology of the level sets of Vh^a-, we 
obtain the characterization of the dynamics of the composition maps on these invariant level 
curves. This result improves the ones obtained in |8] for a6 > 0, where only the stability 
problem was considered. 

We will use the following notations: 



These points correspond to the non-zero fixed points of Gb^a-, that only exist for ab < 1/4. 
We also set /imin := yb,a{P±) = —\ — ab + \fab. 



P± = {x±,y±) : 




Theorem 2. (a) If ah > 1/4, then the origin is the unique fixed point of Gh^a O'l^d it is of 
elliptic type. The sets {Vii^aix,y) = h > 0} surround the origin, are diffeomorphic to 
circles and invariant under G^^a ■ Moreover on each of them the map is conjugated to 
a rotation. 

(h) Ifab < 1/4, then: 

(i) The origin is a hyperbolic saddle of G^^a- The set {V{,^a(x,y) = 0} \ {(0,0)} 
coincides with the stable and unstable manifolds of the saddle and consists of two 
curves diffeomorfic to the real line. Each of them is invariant under G^^a o,nd 
the map is conjugated to a translation. 

(a) The sets {Vb^a{x,y) = h > 0} are diffeomorphic to circles surrounding the set 
{V^b,a < 0}. They are invariant under G^^a CLnd on each of these sets the map is 
conjugated to a rotation. 

(Hi) P± are fixed points of Gb,a of elliptic type. 

(iv) The sets {Vb^aix,y) = h < 0,h ^ /imin} are diffeomorphic to two circles sur- 
rounding and P_ respectively. Each of these circles is invariant under Gb^aj 
and the map on each of these sets is conjugated to a rotation of the circle. 

In next subsection we will use the above result to describe the behavior of the sequences 
{xn\n& generated by the recurrence 

Notice that when the level sets of Vb^a are diffeomorphic to circles, the above result 
allows to associate to each level set {Vb^a{x,y) = h} the rotation number of the map Gb^a 
restricted to it. We denote it by pb,a{h) and we call the function pb^a, the rotation number 
function. Recall that when Pb,a{h) = p/q G Q, with gcd{p,q) = 1, then {Vb^aix,y) = h} 
is full of g-periodic points of Gb^a and when pb,a{h) Q, then the orbit of Gb^a is dense 
on {Vb^a{x.,y) = h}. Notice that pb,a{h) is a well defined function of /i, because when 
{Vb^a{x,y) = h} has two ovals, then pb^a coincides on both because of the invariance of the 
map Gb^a by the change of variables (x,y) — > (— x, —y). Hence, to describe the periods of 
Gb^a we must study the map pb^a{h)- 

Following the techniques introduced in [U [5] and some new tools we prove the following 



Proposition 3. (i) Let Pb,aih) be the rotation number map associated to Gb^a for h > 0. 
Then it is analytic on (0,+oo). Moreover, liuih^-^^ pb^a 

(h) = 1/2 and 



result: 







when ab < 1/4, 




when ab > 1/4. 
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(ii) When ah < 1/4, Pb,a{h) is also defined for h G (/imin) 0) it is analytic on this interval 

and 

lim pf) a{h) = — arccos(l — 16 Vab + 32ab). 

Proof, (i) Following j5l Thm 1] under our hypotheses the rotation number function on each 
level set admits a dynamical interpretation in terms of the time of the flow associated to 
the planar vector field 

Xb,a{x, y) = (^--^' ~^f) ^ ~ '^^y ~ 202; - y + 2xy^). 

More concretely, each closed curve {Vf,^a(a;,y) = /i} is a periodic orbit of period T{h) of the 
differential equation associated to X^^ai and moreover Gb^a[x,y) = ip{T{h);x,y), where (/? 
is the flow associated to Xh^a- Then pb,aih) = T{h) /T{h). From this expression it is not 
difficult to see that it is an analytic function for h € (0, oo), see [^j- 
A straightforward computation shows that when ah > 1/4 



Spec(I?Gb,a(0,0)) = {A+,A-} , where 
By using the tools developed to prove [4, Prop. 19] we get that 



± _ l-2ab±i \/Aab - 1 
~ 2a6 ■ 



lim pb,a{h) = arccos (Re(A)) = -*- arccos ( - 1 
/i-s>0+ ivr ivr \lao 

as we wanted to show. 

To study the case ah < 1/4 we will use another approach. Notice that in this situation 
the level set Vq := {Vi,^a{x, y) = 0}, is formed by the origin, which is a saddle point, and their 
stable and unstable manifolds, which coincide, forming two loops, see Figure 2. Moreover 
on each of these loops the map G^^a is conjugated to a translation and 

lirn G^^„(xi,X2) = (0,0), 

for each (xi,X2) S Vq \ {(0,0)}. Given any e > 0, we take such that 1/(2A^) < e and 
the points := G^^(xi,X2). Observe that both are in the same connected component of 
Vo \ {(0,0)}. Note that Gl^^{P-) = P+ . By continuity of Gf^^ we known that the image 
for G1^ of a point, near P~, on {Vfc^a(x,y) = /i > 0} is, for /i > small enough, as close 
as we want to P^ . This implies that after 2A^ iterates this point has given less that half a 
turn to {Vh^a{x,y) = h > 0}. Consequently, for h small enough, pb,a{h) < 1/{2N) < e, as 
we wanted to prove. 

The case ah = 1/4 follows from the above two results, again by continuity arguments. 
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Finally, to compute lim Pb,a{h) we consider the change of variables given by (z, w) = 
1/y) which conjugates the behavior at infinity of Gf,,a with the behavior at the origin 
of the new map: 

/ z [l + aw'^) p{z,w)w\ 



G{z, w 
where 



l — zw + aw"^^ q{z,w) 



p{z, w) = 1 + 2aw'^ — 2zw + bz'^ + a'^w'^ — 2azw^ + 2abz'^w'^ + z'^w'^ + a'^bz'^w'^, 
q{z, w) = 1 + 2aw'^ - zw + bz"^ + a^w'^ + 2abz'^w'^ + a'^zw^ + aHz^w'^ - az'^w^. 

Clearly G{z,w) = {—z + 02{z,w),—w + 02{z,'w)) and so DG{0,0) has the eigenvalues 
±1 = ibe'^*. Hence arguing as in the study of Gb,a near (0, 0) we obtain that 

lim pb,a{h) = 77- = J, 

as we wanted to prove. 

(ii) Follows by using similar arguments. | 

Remark 4. (i) Applying the method developed in the previous proposition to study Gb we 
would obtain lim Pbih) = 1/4. This gives a new proof of this fact which is essentially 
different of the one given in 0/. 



(ii) It holds that 



1 f 1 \ 1 . ^ 1 

— arccos — - — 1 < — arccos ( — 1) = -. 
2tt \2ab 2-11 ^ ' 2 



Hence the map Pb,a{h) can not be a constant function. In particular this implies that there 
are no positive values a and b for which Gb or Gb,a ctfs globally periodic. 

Before stating our main result on the periods appearing for Gb^a, we study in detail the 
2-periodic points. 

Lemma 5. For a > and b > the map Gb^a has 2-periodic points if and only ifO<ab< 
1/16. Moreover they are given by the two ovals {Vb^aix,y) = —ab}. 

Proof. We have to solve the system Gb,a{Gb,a{x-: v)) = [x, y), which is equivalent to the new 
one Gb,a{x,y) = G^^(x,y). After some manipulations it gives 

J {2x'^y -x + 2by) R{x, y) = 0, 
I (2xy2 + 2ax - y) R{x, y) = 0, 

where R{x, y) := x^y^ + ax^ — xy + by"^ + ab = Vb^a{x, y) + ab. Its solutions have to contain 
also the fixed points of the map. It is not difficult to check that the fixed points are the 
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real solutions of the system {2x^y — x + 2by = 0, 2xy^ + 2ax — y = 0}. Hence the 2-periodic 
points coincide with the set TZ := {R{x,y) = 0}. Therefore we need to study when it has 
real points. Notice R is a quadratic polynomial with respect to x. Its discriminant is 

D{y) := -46/ + (1 - Sab) y^ - 4ba^. 

Since a > and 6 > then TZ will be non-empty when D{y) takes non-negative values. The 
discriminant of D{y), with respect to y^, is 16ab— 1 > 0. So when < a6 < 1/16, the set TZ 
is non empty. When ab = 1/16 it coincides with the other fixed points P±. For ab < 1/16 it 
is formed by two ovals surrounding them, which form precisely the level set with rotation 
number 1/2, that is Pb,a{—o-b) = 1/2. Hence the result follows. | 

Recall that the origin is always a fixed point of the map Gh^a and that when ab < 1/4 
it has two more fixed points. Next result gives the periodic points of this map with period 
greater than one. 

Theorem 6. Let Gh^a the map defined in Q for a > and b > 0. 

(a) For < ab < 1/16 it has continua of periodic points of all the periods greater than one. 

(h) For 1/16 < a6 < 1/4 has continua of periodic points of all the periods greater than 
two and has no periodic points of period two. 

(c) For ab > 1/4 it has continua of periodic points of all the periods q such that there exists 
p GN, with gcd(p, q) = 1 and such that 

- G — arccos — - — 1 , - . 
q \2-K \2ab J 2) 

In particular, there is a computable number qo{a,b) such that Gb^a has continua of 
periodic points of all periods greater than qo{a, b). 

Proof. It is a direct consequence of Theorem [2| Proposition [3] and Lemma [5] Notice that the 
fact that all periods greater than 2 appear when ab < 1/4 is because g > 3, 1/g G (0, 1/2). 
The final conclusion in (c) follows from the well-known fact that any open interval contains 
irreducible fractions with all the denominators bigger that a given number which depends 
on their extremes. A constructive upper bound for this number is developed in | 

Remark 7. (a) If the rotation number function Pb.a{h) was monotonous with respect to h 
then we could also completely characterize all the periods for the case ab > 1/4. They would 
be precisely the ones given in item (c). 

(b) Theorem^ can be easily adapted to obtain similar results for the Gumovski-Mira 
maps, Gb,b > 0. The corresponding results complete the ones of 13] obtained for the case 
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b > 1/2, where all the level sets of the corresponding first integral, but the origin, are 
diffeomorphic to circles. 



2.2 Study of the recurrence ([5j) 

Theorem [2] will allow us to characterize the sequences {x^jneN generated by the recur- 
rence ^ . Observe that in fact since the map Gb^a is a diffeomorphism of , by using its 
inverse we can define the extended sequence {xn}n£Z- Recall that it is said that a recurrence 
is persistent if any sequence {xn}nez is contained in a compact set or M, which depends on 
the initial conditions xi and X2- 

Clearly for any a > 0, b > and xi = X2 = we get that Xn = for all n G Z. Next 
result describes the behaviour of {xn}nez for non-null initial conditions. 

Theorem 8. The recurrence ^ is persistent. Furthermore, let {xn}n& any sequence 
defined by ^ with initial conditions xi and X2, (xi,X2) / (0,0) and set 



-(a^ + b^)-l+ {a + b) ^/ia-b)^ + l 



hmin ■= —- — ab + Vab and /i+ := - 
It holds that hmm </?-+< 0. 

(a) Assume that ab < 1/4. Then: 

(i) For (xi,X2) = P±, which coincide with the set {Vh^a{x,y) = /imin}, the sequence 
{^^njnez is two periodic when a ^ b and a constant sequence when a = b. 

(a) For {xi, X2) ^ (0,0) and Vh a{xi, X2) = 0, lim x„ = 0. 

' n— >itoo 

(Hi) For a numerable and dense set of values ofh, the level sets {Vb^aix,y) = h} give 
initial conditions such that {xn}n&z is a periodic sequence of even period which 
only depends on h. 

(iv) For the rest of values of h (which form a dense set of full measure), the level 
sets {Vb^aix,y) = h} give initial conditions such that the adherence of the cor- 
responding sequence {xn}n£Z is given by one or two intervals. More concretely, 
this adherence is given by one interval when either a = 6; /i > 0; or /i+ < /i < 
and by two intervals when a ^ b and /imin < h < h^. 

(b) Assume that ab > 1/4. Then only the behaviors given in items (Hi) and (iv) appear. 
Moreover in item (iv) the adherence is always formed by a single interval. 

Proof. We will use the techniques introduced in [6] where we studied the 2-periodic Lyness 
recurrences. First notice that the relation between the terms of the recurrence and the 



10 



iterates of the composition map is given by 



Gb,a{x2n-l,X2n) — ix2n+l, X2n+2) , Ga,bix2n, X2n+l) — {x2n+2, X2n+3) , 

where {xi,X2) € and n G Z. Also observe that a simple computation shows that 

Vb,a{x,y) = Va,b{Ga{x,y)), 

which imphes that the odd terms of the sequence {xn}nez are contained in the projection 
on the X-axis of the level sets of the form {Vb^a{x, y) = Vb^aixi, X2) = h} and the even ones 
are in the corresponding projection of the level sets {Vafi{x,y) = Va^biGa{xi, X2)) = h}. 

The result (i) is clear. Statement (ii) is a new consequence of Theorem [2] and the fact 
that the sets {Vb^a{x,y) = 0} \ {(0,0)} are two homoclinic loops starting and ending at 
(0,0). To prove (iii) and (iv) we consider the cases (a) and (b) together. 

Again from Theorem [2] we know that the sequence starting at a point (xi, X2) such that 
Pb,aiyb,a{xi,X2)) = p/q G Q with gcd(p, (?) = 1, will be 2g-periodic. When it holds that 
Pb,aiVb,aixi, X2)) Q, then it will densely fill the projection of {Vb^a{x,y) = h} (the odd 
terms) together with the projection of {Va,bix,y) = h} (the even terms). 

Hence to end the proof we only need to show that the rational values of Pa,b(yb,aih)) are 
taken only for a numerable dense set of values h and, when Pa,biy^b,aih)) Q, to distinguish 
whether the above to projections do or do not overlap. 

Let us prove first that /5a,fe(^b,a(^) G Q only for a numerable and dense set of values 
of h. First recall that from Proposition |3] and Remark |4] we know that the function is 
analytic and nonconstant in the intervals where it is defined. Next we use the well-known 
fact that a nonconstant analytic function has a numerable number of zeroes. This implies 
that for any r € Q the set 7^,. := {/i G M : Pb,a{h) — r = 0} is numerable. Therefore Ur-eQ'Pr 
is numerable as we wanted to see. That it is dense follows from the continuity of all the 
involved functions. 

To end the proof we only need to decide whether the projections into the x-axis of 
{yb,a{x-,y) = h} and {Vafi{x^y) = h} intersect or not. 

When /i > 0, for any positive a and h both ovals {Vb^a{x,y) = K} and {Va^b{x-,y) = h} 
surround the origin. Hence the adherence of {xn}ri& is always one closed interval. 

Assume now that /i < (and so ah < 1/4). An easy computation shows that when 
a ^ b, Ga{P±) is neither P± nor Pzp, so the points P± give rise to two different 2-periodic 
solutions of the recurrence. When a = b, P± are fixed points of Ga, and the situation is like 
the previous one. So in the following we also assume that a ^ b. 

Since Ga{P±) is neither P± nor Pz^, and Ga is a diffeomorphism, it is possible to take 
ovals of the form {Vb^a{x,y) = h} with hmin ^ h < whose projections into the abscissa 
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axis do not overlap with the corresponding ones of {Vafi{x, y) = h}. However, by increasing 
h the projections wih start to overlap. 

Let be = I^(a,b, h), and = J^{a,b,h) be the projections on the abscissa axis 
of the ovals contained in the level sets given by {Vb^aix,y) = h} and {Va^b{x,y) = h} 
respectively, surrounding P+ and P_. We want to detect the values of h for which the 
intervals and J"*" (and /~ and J~ respectively) have exactly one common point. First 
we seek for their boundaries. Since the level sets are given by quartic curves, quadratic 
with respect the y-variable, these points will correspond with values of x for which the 
discriminant of the quadratic equation with respect to y is zero. So, we compute 

Ri{x,h,a,b) := dis {Vb,aix,y) - h,y) = -4ax^ + {-iab + ih + 1) .x^ +4bh, 
R2{x,h,a,b) := dis (V;,b(x, y) - h,y) = -Abx"^ + {-4ab + Ah + 1) x'^ + 4ah. 

To search for relations among a,b and h for which these functions have some common 
solution, X, we compute the following resultant 

Res{Ri,R2;x) = 256 (a - b)^ h'^ (l6h^ + (iGa^ + 16 6^ + 8) /i + (46a - 1)^)^ 

Since a ^ b, the possible bifurcation values of h are given by the roots of the quadratic 
polynomial 

Pi{h) :=16/^^ + (l6a^ + 166^ + 8)^ + (46a-l)^ (7) 



which are 



-(a' + ft') - ^ ± (a + 6) Via-b^ + l 



Then the number of disjoint intervals depends only on the relative positions of /iminj h± and 

Vb,a{xi,X2). 

To study the relative positions of /imin and h±, we introduce the auxiliary polynomial 
P2{h) = {h + ab + j)^ — ab, because one of its roots is /imin- Then 

Res(Pi, P2; h) = 16 {Aab - if (a - bf. 

This means that we can split the parameter set A := {(a, 6), a,b > 0, ab < 1/4} into the 
following disjoints sets: 

A= {a,b > 0, ab < 1/4, a > bjU {a,b > 0, ab < 1/4, a = b} U {a, 6 > 0, ab < 1/4, a < b}. 

Observe that, because of the properties of the resultant, the relative positions of /imin and 
h± on the above subsets do not vary. So testing their values for some particular parameters, 
we obtain that when a / 5 we always have that /imin < /i+ < (we also notice that when 
a = b, h- < /imin = h+). 
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In summary, the above considerations imply that when a ^b, and since Ga{P±) 7^ P±, 
the projections into the abscissa axis of all the ovals {Vf,^a(x,y) = h} and {Va^b = h} when 
^min < /i < /i+ do not Overlap, giving rise to two different intervals. When /i+ < /i < they 
overlap giving only one interval. | 

From the above theorem we get that the list of behaviours (i)-(iv) listed in the introduc- 
tions are the only ones that the sequences generated by ([5]) can present. Moreover, given 
some values of a and b the information about the possible periods of the periodic sequences 
can be obtained from Theorem |6l 

We end this subsection with a concrete example. 

Example 9. Consider a = 0.01, b = 0.49. For any initial condition (xi,X2), set Hq = 
yb,a{xi,X2)- From the definitions given in Theorem^ we have that /imm = —0.1849 and 
h+ ~ -0.0928, so 

• If hQ G (/iinm)^+)) ihe sequence generated by the recurrence ^ is either periodic (with 
even period) or it fills densely two intervals. 

• If ho €z [h-f-,0), the sequence generated by the recurrence ^ is either periodic (with 
period a multiple of 2) or it fills densely one interval. 

For instance if we consider the initial conditions (xi, X2) = (x^ — 0.1, y+ — 0.1), we have 
that ho = Hb^a{xi, X2) — —0.1459, hence ho G {hniin,h+), and the corresponding sequence 
generated by the recurrence ^ is either periodic or it densely fills two disjoint intervals. 
This last option is the one that appears in our numerical simulations, see Figure 3. 























If 


























. i 















Figure 3: Here a = 0.01, b = 0.49. The ovals {Vb,a{x,y) = ho} and {^4,6 = ^0} with 
ho — —0.1459 and their projections on the abscissa axis. 
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2.3 Integrability for general /c-periodic Gumovski-Mira recurrences 

In this section we search for invariants of the form 

4 

I(x,y,n)= Ilfx'yi (8) 
i+j=i 

for the recurrence ([5]), but for any sequence {(in}n of positive numbers. This form is 
natural because includes the invariants known for {/3n}n constant or 2-periodic. As far 
as we know, this method was introduced in [T^. It was also used in |7j to prove that 
non-autonomous Lyness recurrences, Xn+2 = («ra + Xn+i)/xn-, only admit invariants of a 
suitable form when {an}n are periodic with periods 1,2,3 and 6. If we could find a new 
integrable Gumovski-Mira recurrence all the methods used in this section could be adapted 
to study it. Unfortunately, next result shows that for the Gumovsky-Mira non-autonomous 
recurrences, only the cases of {/3ra}n with period 1 and 2 seem to be integrable. We have 
chosen /(x, y, n) with terms of degree at most 4 for simplicity. Bigger degrees could be also 
studied by using the same approach. 

Proposition 10. The non- autonomous k-periodic recurrence with a sequence of 

positive numbers, has invariants of the form ^ if and only if {f3n}n is k-periodic and k G 
{1, 2}. Moreover, in these cases the invariants are I{x, y, n) = jinx"^ + Pn+ilP +1x^11^ — ixy , 
for any 7 7^ 0. 

Proof. The condition that a function I{x, y, n) of the form ([S]) is a non-autonomous invariant 
of the recurrence ([s]) writes as 

l[y,-x+ ^ n + 1 ) - I{x, y, n) = 0, 

for all {x, y) G and all n G N. Imposing that all the monomials of the left hand side of 
the above relation vanish and after some work we get in particular that /3n+2 = Pn and so 
that such an invariant exists only when {/3n}n is constant or 2-periodic. For both cases we 
also obtain the explicit form of the known invariants. | 



3 The chaotic case 



In this section we consider the non-autonomous recurrence 



Xn+2 = -Xn + — n , (9) 

l + xt 



with 

a > for n = 2k, 



OLr. 



b>0 for n = 2k + l, G N. 
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Recall that for each fixed value of a, the map associated to the corresponding au- 
tonomous recurrence 

Xn+2 = -Xn + 2 ' (10) 



has the first integral 



Wa{x,y) = xy +x +y - axy. 



Moreover it is conjugated to Gi/a- Therefore Fa and Fa^a are integrable and the dynamical 



systems generated by them as well as the recurrences generated by (10) can be easily 
described using the results of the previous section. 

3.1 The map as a perturbed twist 

Here we will see that the map -Fft^a; which writes as, 

r f f , b{l + y^){ay-x{l + y^)) \ 
Fb^a{x,y) = [-X + — — 2^~y+~n~, — 2^2 , / rT, — 2v^ ' 

is far from preserving neither the integrability nor the translation- like dynamics, proved for 
Gb^a- In fact, our numerical simulations show all the features of a non-integrable perturbed 
twist map, that is: many invariant curves and, between them, couples of orbits of p-periodic 
points (for several values of p) , half of them of elliptic type and the other half of hyperbolic 
saddle type, see for instance [21 Chapter 6]. Notice that the sequences corresponding to 
the stable manifolds of the hyperbolic saddles are precisely the ones that give rise to the 
behaviour (v) described in the introduction. 

More precisely, for instance, when a £ (0, 2) the positive level sets of Wa in are 
closed curves surrounding the origin, which is the unique fixed point of Fa (see Figure 1 
(a)) and we know that the action of Fa restricted to each of them is conjugated to a rotation. 
In fact, both Fa and Fa can be seen as area preserving integrable twist diffeomorphisms 
since, at least locally, the rotation number function associated Fa to most level set of Wa is 
monotonous. Therefore when we consider Fa+e,a, for e small enough, it can be thought as 
a perturbed twist map. Our numerical experiments show that this perturbed twist map is 
no more integrable and that this behaviour remains even when the values of b are not close 
to the ones of a. This is the situation depicted in Figure 1 (b) and (c). Unfortunately we 
have not been able to prove the non-integrability of F^^a for b ^ a. 

When a > 2 the negative level sets of Wa are two nests of closed curves surrounding two 
fixed points of Fa- These two nests are surrounded by a poly cycle composed by a saddle 
point of the map and a couple of homoclinic curves surrounding both elliptic fixed points. 
The rest of the level sets are given by closed curves surrounding the polycycle, see Figure 2. 
The dynamics of Fa on the invariant level sets of Wa is again translation-like. In this case. 
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when we consider b = a + e for e small enough, near the nests of Fa the map -F^^a can 
be seen again as a non-integrable perturbation of the twist map = Fa^a, and the same 
phenomena than above appear. These behaviours are also observable when the values of b 
are not close to the ones of a. This is the situation exemplified in Figure 4, where it can 
also be observed orbits of -F^ ^ which seem to fill densely a sets of positive Lebesgue measure 
in M^. 




(a) (b) 

Figure 4: (a) Several orbits of the map ^9/2,4- (b) Plot of 10^ iterates of the orbit of Fg/2,A 

starting at (1.18,0.1). 

From the viewpoint of the solutions {xn}n& of the equation Q, another difference is 
that the number of disjoint intervals of the adherence of an orbit can be more than 2, as it 
is shown in the following examples. 



Example 11. Consider a = 2 and b = 1/2. We give some initial conditions giving rise to 
solutions {xn}nez which, numerically, seems to densely fill more than 2 intervals. 



Initial conditions 


Intervals 


(1.10,0.5) 


7 


(1.25,0.5) 


16 


(1.29,0.5) 


1 


(1.30,0.5) 


23 


(1.35,0.5) 


1 


(1.40,0.5) 


6 


(1.48,0.5) 


6 



16 



In Figure 5 are depicted the 6 intervals corresponding to the 
and the 16 intervals of the initial condition (1.25,0.5) 



initial condition ( 



1.48,0.5), 
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(a) 



Figure 5: Orbits of i^i/2,2 ^-^id their images under F2 corresponding to the initial 
conditions (1.48,0.5) and (1.25,0.5) respectively, together with their corresponding 

projections on the abscissa axis. 



□ 



Example 12. An example similar to Example 11 with a = 4 and h = 9/2. 



Initial conditions 


Intervals 


(2.5,0.5) 


4 


(2.9,0.5) 


1 


(3.0,0.5) 


5 


(3.1,0.5) 


3 


(3.2,0.5) 


1 


(3.6,0.5) 


4 



In Figure 6 we show the 4 intervals corresponding to the initial condition (2.5,0.5), and the 
5 intervals of the initial condition (3,0.5). 
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(b) 



Figure 6: Orbits of -^9/2,4 and their images under F4 corresponding to the initial conditions 
(2.5,0.5) and (3,0.5) respectively, and their corresponding projections on the abscissa axis. 



3.2 Local study of Fb^a near the origin 

The origin is an elliptic fixed point of the area preserving map F^^a when a6 < 4. In fact 
the characteristic polynomial of DF^^ai^^O) is 

Pa,bW = + (2 - ba)X + 1 



and Spec(i:>Fb,„(0,0)) = {A+, A"}, where A± := -1 + f ± ^ y^ab{4.-ab). 
Computing the resultants 

Res(Pa,;,(A),A'' - 1; A), for /c = l,2,3 

we know that for ab 1 the eigenvalues A^ are under the non-resonance conditions of 
Moser's Stability Theorem, see [H [2D]. That is, they are not /c^^-roots of the unity for 
k = 1,2,3. We remark that Example 11 seems to show that the situation a6 = 1 is not 
special at all. 

We have obtained the local Birkhoff normal form of -F^ ^ at the origin using [12\ [T9] . It 
is given by 2: — )■ X'^z (1 + iazz + 03(2;, z)) , where z = x + iy G C and 



_ 3 (6 + o) y/gb (4 - ab) 
~ 46 (4 - ab) ^ 

So, from Moser's Theorem we know that for ab 7^ 1, in every neighborhood of the 
origin, there exist infinitely invariant closed curves surrounding it. Although the theorem 
does not allow to make the distinction between the integrable and the non-integrable cases, 



18 



by using the Poincare-BirkhofF Theorem, see |20j, we can ensure the existence of periodic 
orbits having ah periods bigger that a given one, between each two of these invariant 
closed curves. Therefore we have proved that for the corresponding values of a and b the 
recurrence ([9| generates periodic sequences of any even period, bigger that some constant, 
which depends on a and b. These sequences correspond to the ones described in item (i) of 
the introduction. 
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